Prepared for submission to JHEP 



A Reference Model 
o : for Anomalously Interacting Bosons 

a 



> 

00 



in 

o 
o 



in 

Dh; M. v. Chizhov'^'' 

^~H, "^CSRT, Faculty of Physics, Sofia University, 1164 Sofia, Bulgaria 

r^ I DLNP, Joint Institute for Nuclear Research, 141980, Duhna, Russia 



E-mail: mih@phys.uni-sofia.bg 

Abstract: A simple reference model for anomalously interacting bosons is proposed and 
Ps| . implemented in the CompHEP package. This allows preparing for an experimental search 

of these bosons at powerful colliders, such as Tevatron and LHC. New signatures and some 
experimental consequences are shortly considered. 

Keywords: excited bosons: Zstar and Wstar 



X ■ ArXiv ePrint: 1005.4287 

M'. 

Contents 

1 Introduction 2 

2 The reference model and its implementation 2 

3 Consequences for colliders 4 



1 Introduction 

The hypothetical heavy neutral bosons, Z', interacting with the known fermions, quark 
and leptons, should have a clear signature at the hadron colliders. Due to their resonance 
production a peak in the dilepton invariant high-mass distribution of the Drell-Yan process 
can be expected. In the case of leptophobic bosons their decay in the heavy quarks, t and 
b, allowing tagging and almost a full kinematics reconstruction, can be investigated. 

The extra gauge bosons, Z', occur in any extensions of the Standard Model (SM) gauge 
group. If such extensions are motivated by the Hierarchy Problem, they lead to a new 
physics around TeV energies not far above the weak scale. Needless to say, understanding 
experimental consequences of the latter is of a fundamental importance. 

However, usually only the gauge bosons, which have minimal couplings with the SM 
fermions 

^z' = E (sL i^ir^L + aL i^iin^ii) ^ , (i.i) 

/ 

are considered in literature. In the recent paper [1] it was shown, that there are at least 
three different classes of theories, all motivated by the hierarchy problem, which predict also 
the appearance of new spin-1 weak-doublets, V^ = {Z^,W*), with the internal quantum 
numbers identical to the SM Higgs doublet. 

In contrast to the well-known Z' and W bosons they possess only effective anomalous 
(magnetic moment type) couplings with ordinary light fermions [2] 

^y = JlY. (^L^^'^^^ ^[.K] + gL Dl Vl^ Vna^'^i) , (1.2) 

/ 

where M is the scale of the new physics, "0^ are the right-handed singlets and \I'^ are the 
left-handed doublets. Z?^ are the usual SU{2)w x C/(l)y covariant derivatives, and the 
obvious group and family indexes are suppressed. 

Although many Beyond the Standard Models (BSM) have been implemented in Monte 
Carlo (MC) Generators like PYTHIA, the new bosons are not yet even assigned to any 
number within the MC particle numbering scheme. The aim of this paper is to suggest a 
simple reference model and implement it in one of the most reliable among Matrix Element 
MC Generators, the CompHEP package [3]. This will allow to perform simple calculations 
of the processes involving the new bosons. In agreement with Les Houches Accord [4] the 
corresponding generated event files can be passed to PYTHIA for showering/hadronization 
and further also for simulation and reconstruction within the real experimental software 
framework. 

2 The reference model and its implementation 

In order to construct the simple reference model for experimental verifications we need to 
derive the main features of the new spin-1 bosons from a whole variety of possible models. 
For example, the doublet structure always introduces new particles in pairs: neutral Z* 
and charged W* spin-1 bosons. Moreover, due to non-zero hypercharge of the doublet the 



Table 1. Independent parameters of the _ESM model (varsX.mdl) 



Variables 



Name 


1 Value 


> Comment 


<l 


GW 


0.652 


new coupling constant = EE/SW 




MZd 


1000. 


mass of down-type Z* 




MZu 


|700. 


mass of up-type Z* 




MWX 


500. 


mass of W* 




sx 


|0.5 


sin of charged boson mixing 





Table 2. List of particles of the _ESM model and their properties (prtclsX.mdl) 



Particles 



Full name 


IP 


|aP 


|2* 


spin 


mass 


1 width 


color 


aux 


1 > LaTcX(A) 


< 


> LaTcX(A+) 


<l 


down-type Z* 


|Zd 


Zd 


2 




|MZd 


|wZd 


|1 




Z-*_d 




|Z"*_d 




up-type ReZ* 


|Zr 


|Zr 


2 




|MZu 


|wZr 


|1 




ReZ"*_u 




|ReZ"*_u 




up-type ImZ* 


Zi 


|Zi 


|2 




|MZu 


|wZi 


1 




lmZ"*_u 




|lmZ"*_u 




W* 


X+ 


|x- 


|2 




|MWX 


|wWX 


1 




|w-{*+} 




|w-{*-} 





minimal particle set consists of the four spin-1 bosons: the two charged W states and 
the two neutral CP-even ReZ* and CP-odd ImZ* states. 

The crucial common feature of all approaches considered in [1] is, that the lightness 
of the Higgs doublets is guaranteed, because they are related to the spin-1 V^ bosons by 
symmetry. It means that each spin-0 Higgs doublet is associated with the corresponding 
spin-1 doublet and vice versa. In the supersymmetric SM extensions at least two Higgs 
doublets are needed, which demands introduction of more than one new spin-1 weak- 
doublet. In order to prevent the appearance of flavor-changing neutral currents one doublet 
should couple only to up-type quarks while the other should couple to down-type quarks 
and charged leptons only [5]. 

The different doublets are associated with the different vacuum expectation values of 
the Higgs doublets and, in general, acquire different masses. Let's choose them according 
to the estimations in [2] 



Mz* ~ 700 GeV, 



M^. ~ 



1000 GeV. 



(2.1) 



In contrast to the case of neutral bosons the charged bosons can mix, which results in 
an additional mass splitting. So, the mass of the lightest boson can be estimated as 
Mw* ~ 500 GeV [2] . It will be shown that all properties of the new bosons can be described 
by the following independent input parameters: coupling constant, masses and mixing of 
the charged bosons. All input parameters of the model called _ESM are summarized in 
Table 1. 

Table 2 includes all characteristic representatives of the _ESM model and their prop- 
erties. 



Table 3. Parameters depending on the basic ones (funcX.nidl) 



Constraints 



Name 


1 > Expression 








< 1 > Comment < | 


CX 


sqrt(l-SX-2) 








|cos of charged boson mixing 


pi 


|acos(-l) 








|Pi 


Sqrt3 


|sqrt(3) 








|sqrt of 3 


rtZ 


Mtop-2/MZu-2 








|(Mtop/MZu)-2 ratio 


rtW 


|Mtop-2/MWX~2 








|(Mtop/MWX)-2 ratio 


wZd 


|GW~2*MZd/(4*pi) 








|width of down-type Z* 


wZr 


GW~2*MZu/(12*pi) 


*(2+(l 


+8*rtZ) *sqrt (l-4*rtZ) ) 


|width of up-type ReZ* 


wZi 


|GW-2*MZu/(12*pi) 


*(2+sqrt((l- 


-4*rtZ)-3)) 


|width of up-type ImZ* 


wWX 


|GW-2*MWX/(4*pi)* 


(1+sx- 


2*rtW 


■2*(3-2*rtW)/12) 


1 width of W* 



Their interactions can be derived from (1.2) 



-ref 



^ {da^'^d + ecj^^e) df^Z^i^ + /^^ {ua^'% d^ReZ*^^ + i ua^'^j^u df,lmZ*^^) 



V2Mz:; 



V^Mz 



+jf^ Un9x Wcj'"' dR + -^ cos OxUR^^'dL + sin 6 xiTZa^^'eR ) a^Wl+ + h.c, (2.2) 



where for siinphcity the scale of the new physics is equal to the mass of the corresponding 
boson and g being the coupling constant of the SU{2)\y weak gauge group. The flavor 
universality is assumed and the family indexes are omitted. The coupling constants are 
chosen in such a way that all fermionic decay widths of the new bosons are the same in the 
Born approximation and in the massless quarks limit. This value coincides with the total 
decay width of the heavy W' boson of the same mass. The corresponding decay widths are 
presented in Table 3 taking into account the top quark mass only. 

The present paper discusses only the resonance production of the heavy bosons and 
their subsequent decay into fermion pairs, which is important for early discoveries. In 
the case of the light fermions it is impossible to discriminate the multiplicative quantum 
numbers of the neutral boson, P and C, because they have identical signatures. Therefore, 
we will consider only one of the down-type neutral bosons, for instance, HeZ^. Table 4 
describes the vertex structures of the new interactions (2.2). 

These four tables completely describe the model and can be used for process evaluations 
within the CompHEP package. 



3 Consequences for colliders 



The hadron colliders, due to their biggest center-of-momentum (CM) energy t/s ~ several 
TeVs and their relatively compact sizes, still remain the main tools for discoveries of very 
heavy particles. Since the production mechanism for new heavy bosons at a hadron collider 
is the qq resonance fusion, the presence of partons with a broad range of different momenta 
allows to flush the entire energetically accessible region, roughly, up to \/s/2. This, in a 
way, fixes the dominant production and decay mechanisms. 



Table 4. Vertices of interactions of the _ESM model (IgrngX.mdl) 



Lagrangian 
PI 
B 

D 
E 

L 

n 
s 
c 

T 
U 
C 
T 
U 
B 
C 
D 
E 
L 

n 

Ne 

Nl 

Nm 

S 

T 

U 



|P2 


|P3 


P4 1 


|b 


|Zd 




|d 


|Zd 




e 


|Zd 




1 


|Zd 




m 


|Zd 




s 


|Zd 




c 


|Zr 




t 


|Zr 




u 


|Zr 




c 


|Zi 




t 


|Zi 




u 


|Zi 




t 


|x- 




s 


|X+ 




u 


|x- 




nu 


|x- 




Inl 


|x- 




im 


|x- 




e 


|x+ 




1 


|x+ 




m 


|x+ 




c 


|x- 




|b 


|X+ 




|d 


|X+ 





I > Factor < 

|GW/(2*Sqrt2*MZd) 

|GW/(2*Sqrt2*MZd) 

|GW/(2*Sqrt2*MZd) 

|GW/(2*Sqrt2*MZd) 

|GW/(2*Sqrt2*MZd) 

|GW/(2*Sqrt2*MZd) 

|GW/ (Sqrt2*Sqrt3*MZu) 

|Gtf/(Sqrt2*Sqrt3*MZu) 

|GW/(Sqrt2*Sqrt3*MZu) 

|GW/ (Sqrt2*Sqrt3*MZu) 

|GW/ (Sqrt2*Sqrt3*MZu) 

|GW/ (Sqrt2*Sqrt3*MZu) 

|GW/(4*Sqrt3*MWX) 

|GW/(4*Sqrt3*MWX) 

|Gtf/(4*Sqrt3*MWX) 

|GW/(4*MWX) 

|GW/(4*MWX) 

|GW/(4*MWX) 

|GW/(4*MWX) 

|GW/(4*MWX) 

|GW/(4*MWX) 

|GW/(4*Sqrt3*MWX) 

|GW/(4*Sqrt3*MWX) 

|GW/(4*Sqrt3*MWX) 



> dLagrangian/ dA(pl) dA(p2) dA(p3) 
|G(m3)*G(p3)-G(p3)*G(m3) 
|G(m3)*G(p3)-G(p3)*G(m3) 
|G(m3)*G(p3)-G(p3)*G(m3) 
|G(m3)*G(p3)-G(p3)*G(m3) 
|G(m3)*G(p3)-G(p3)*G(m3) 
|G(m3)*G(p3)-G(p3)*G(m3) 
|G(m3)*G(p3)-G(p3)*G(m3) 
|G(m3)*G(p3)-G(p3)*G(m3) 
|G(m3)*G(p3)-G(p3)*G(m3) 

(G (m3) *G (p3) -G (p3) *G (m3) 
I (G(m3) *G (p3) -G(p3) *G (m3) 
I (G(in3) *G (p3) -G(p3) *G (m3) 

(G(in3)*G(p3)-G(p3)*G(in3) 

(G(m3)*G(p3)-G(p3)*G(m3) 

(G(in3)*G(p3)-G(p3)*G(in3) 
|(G(in3)*G(p3)-G(p3)*G(iii3) 
I (G (m3) *G (p3) -G (p3) *G (m3) 
I (G(m3) *G (p3) -G(p3) *G (m3) 
|(G(m3)*G(p3)-G(p3)*G(m3) 
|(G(m3)*G(p3)-G(p3)*G(in3) 
|(G(m3)*G(p3)-G(p3)*G(m3) 

(G(iii3)*G(p3)-G(p3)*G(m3) 
I (G(m3) *G (p3) -G(p3) *G (m3) 
I (G(in3) *G (p3) -G(p3) *G (m3) 



*G5 

*G5 

*G5 

*(Sqrt3*SX*(l-G5)+2*CX*(l+G5)) 

*(Sqrt3*SX*(l+G5)+2*CX*(l-G5)) 

* (Sqrt3*SX* (1-G5) +2*CX* ( 1+G5) ) 

*SX*(1-G5) 

*SX*(1-G5) 

*SX*(1-G5) 

*SX*(1+G5) 

*SX*(1+G5) 

*SX*(1+G5) 

*(Sqrt3*SX*(l-G5)+2*CX*(l+G5)) 

*(Sqrt3*SX*(l+G5)+2*CX*(l-G5)) 

*(Sqrt3*SX*(l+G5)+2*CX*(l-G5)) 



In paper [6] it has been found that tensor interactions lead to a new angular distribution 
of the outgoing fermions 



da{qq ^ Z*/W* ^ //) 



dcos0 
in comparison with the well-known vector interactions 

da{qq^Z'/W'^ff) 



oc cos 0, 



dcos0 



oc 1 + cos 



(3.1) 



(3.2) 



It was realized later [2] that this property ensures distinctive signature for the detection of 

the new interactions at the hadron colliders. At a first glance, the small difference between 

the distributions (3.1) and (3.2) seems unimportant. However, the absence of the constant 

term in the first case results in new experimental signatures. 

First of all, since the angular distribution for vector interactions (3.2) includes a 

nonzero constant term, this leads to the kinematical singularity in pT distribution of the 

final fermion 

1 

oc 



cos 9 

in the narrow width approximation T <^ M 

1 



{M/2f 



Pt 



(3.3) 



vr 



(s - M2)2 + Af 2^2 MT 



6(8 -M^). 



(3.4) 



This singularity is transformed into the well-known Jacobian peak due to a finite width of 
the resonance. In contrast to this, the pole in the decay distribution of the Z* /W* bosons is 
canceled out and the fermion px distribution reaches even zero at the kinematical endpoint 
Pt = M/2. Therefore, the Z* /W* boson decay distribution has a broad smooth hump 
with a maximum below the kinematical endpoint, instead of a sharp Jacobian peak (the 
left-hand plot of Fig. 1).^ 




100 200 300 400 500 600 700 800 

P, [GeV] 




Figure 1. Left: the final electron px distributions from the Z^ (solid) and Z'ggj^j (dotted) bosons 
decays. Right: the differential cross-sections for the Z* (solid) and Z'ggj^j (dotted) bosons as 
functions of the pseudorapidity of the final electron. 



According to Eq. (3.f), for the Z* /W* bosons there exists a characteristic plane, per- 
pendicular to the beam axis in the parton rest frame, where the emission of the final-state 
pairs is forbidden. The nonzero probability in the perpendicular direction in the labora- 
tory frame is only due to the longitudinal boosts of the colliding partons. This property is 
responsible for the additional dips in the middle of the final fermion pseudorapidity distri- 
butions of the anomalously interacting bosons, in contrast to the Z'/W' bosons with the 
minimal gauge couplings (the right-hand plot of Fig. f ). 

While the maximum of the gauge boson distribution is centered at the small lepton 
pseudorapidities, which correspond to the central part of the detector, the chiral boson 
distribution has minimum in this region and its maxima are placed at the edges of the 
CDF and DO central calorimeters \r]\ ~ 1. Based on the fact, that the major part of the 
leptons, stemming from the Z' decays, are emitted in the central detector region, both 
collaborations have analyzed the spectrum of the transverse high-energy electrons only in 
the central electromagnetic calorimeters |??e| ^ Vcut — 1- The acceptances for the processes 
with the chiral Z^ bosons and the sequential Z'ggj^j bosons are shown in the left-hand plot 
of Fig. 2. 

The cuts in the backward-forward regions lead to a miss in an essential part of the 
events from the chiral boson decays due to the mentioned previously specific angular distri- 



^ All calculations are performed for the resonant production of the down- type chiral Z*i bosons with the 
mass M = 890 GeV and for the sequential Z'ggi^j bosons with the mass M — 1000 GeV and their subsequent 
decays into e'^e~ pairs in pp collisions at ^/s = 1.96 TeV using the CompHEP package 4.5.1. Thus the 
chosen masses lead to the identical cross sections of 5.6 fb. 



a 0,8 



^ 0,4 
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Figure 2. Left: the acceptances for the Z^ (solid) and Z'ggj^j (dotted) bosons as function of 
\Ve\ < 'ycut cut. Right: the ahowed region (white) for the couphng constant— mass relation of the 
Z^ resonance from the DO Tevatron constraints [7] . The horizontal line corresponds to the ordinary 
value of the SU{2)w coupling constant g. 

bution. As seen from the left-hand plot of Fig. 2 the curve for the chiral Z^ boson mostly 
lies under the Z'^gj^j- curve, and at r]cut — 1 there are around 65% detected events in the 
case of the sequential Z'^gj^j bosons decays and only 45% in the case of the chiral Z^ bosons 
decays. 

Taking into account this fact we can correct and apply experimental constraints of the 
DO Collaboration (TABLE II from [7]) on production cross section in the case of the chiral 
Z^ bosons with arbitrary coupling constant g (the right-hand plot of Fig. 2). One can see 
that due to a smaller acceptance for the Z^ bosons than for the Z'g^j^j bosons, even for the 
same production cross sections more lighter Z^ masses are allowed. For example, in the 
case when the coupling constant g is equal to the ordinary value of the SU{2)\y coupling 
constant, we obtain Mz* > 940 GeV, in comparison with DO result M^' > 1023 GeV. 
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